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Abstract 

We give a general expression for the static potential energy of the gravitational inter- 
action of two massive particles, in terms of an invariant vacuum expectation value of the 
quantized gravitational field. This formula holds for functional integral formulations of eu- 
clidean quantum gravity, regularized to avoid conformal instability. It could be regarded as 
the analogue of the Wilson loop of gauge theories and allows in principle, through numerical 
lattice simulations or other approximation techniques, non perturbative evaluations of the 
potential or of the effective coupling constant. 



1 Introduction. 

The present paper is concerned with the problem of the energy of the gravitational field. This energy has 
been under investigation since the birth of General Relativity and some issues, like the determination of the 
total energy of a field configuration, have been settled in a rigorous way in the ADM formalism or through 
Noether's theorem [||. Other points, however, like the possibility of "localizing" the gravitational energy, are 
still obscure (see [|j and references therein). 

*A. Von Humboldt Fellow. 



Here we shall concentrate on a special case, namely that of the static potential energy of the gravitational 
interaction of two massive bodies. The gauge-invariant description of this interaction in terms of quantum fields 
is a crucial issue for any physically reasonable theory of quantum gravity. 

This potential energy will turn out to be related to the vacuum average of a simple gauge- invariant functional 
of the field, a kind of correlation between "scalar" Wilson lines. While evaluation of this average for a weak field 
on a flat background yields the Newton potential in the usual form , non-perturbative evaluations could show 
modifications in the effective coupling constant or in the dependence of the potential on the distance between 
the particles. 

Like in the case of the Wilson loop, our formula can be implemented quite naturally on a lattice version 
of the theory, in order to allow numerical simulations (Hamber and Williams, |4|). The first results are very 
interesting, showing that in the strong-gravity region of the phase space of lattice theory the potential is indeed 
Yukawa-like, with a mass parameter which decreases towards the critical point (where the average curvature 
vanishes). The effective value of Newton constant has been estimated too. 

This paper extends and refines an earlier work ||]. The outline is the following. In Section 2 we recall 
in brief the general-relativistic view about the gravitational potential energy. In Section 3 we give a formula 
which, treating two masses as external sources for a quantized gravitational field on a flat background, allows 
to write the potential energy of their interaction. This is done using a known technique of euclidean quantum 
field theory ^ . In Section 4 we generalize the formula to the case of "strong" gravity, introducing a definition 
of the source which avoids the use of background metric. Finally, in Section 5 we discuss the compatibility of 
the background formula with the general one and the possible behaviour of the effective coupling constant. 



2 Potential Energy Versus ADM Energy. 

In this Section we work in (3+1) dimensions and follow the conventions of Weinberg |Q. Starting from next 
Section, the metric will be euclidean. The Einstein action is given by 
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S'Einst. = - J^:^ I d'^x ^/g{x) R{x) (1) 



and the action of a material particle of mass m is 



5'Mat. ^ -m dp ^ -g^i,[x{p)]xt^{p)x'' {p), (2) 



where x^{p) is the trajectory of the particle and p is any parameter. The dots will always denote differentiation 
with respect to the parameter. Finally, the Einstein equations have the form 

Rfj^f - ^9t^iyR = -SttGT^i,. (3) 

Let us now decompose the metric in the traditional way 

9^L,^{x) = rjf^i^ + hf_,^{x) (4) 

and denote the linearized Einstein equations in the harmonic gauge as 

-^fZ ^pa — T^i,. (5) 

The inverse of the kinetic operator K is the well-known Feynman-De Witt propagator 

^n^pAx-y)^ —— / „ ,^2 ■ ■ (6) 
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Let us compute the Newton potential starting from the preceding equations. This is an elementary cal- 
culation, but the integral one encounters will be useful in the following. The field produced by a generic 
four-momentum source Tp„ in the linearized approximation is given by 

Ku{x) - j d^v [K-Xu{x - y) TpAy); (7) 

when the source is a particle of mass m at rest in the origin, the only non vanishing component of Tp^ is 

Too{y)^md^y), (8) 

so we have 

liOO/ 



= m r°°dy° [ d^y ^^'(y) 



J-00 ^ - (a;" - y^r - 1^1 

This is the correct result, since in the newtonian approximation we have 

(?oo = -l-2F, (10) 

where V is the Newton potential. 

More generally, one can use the ADM formula. In classical General Relativity the total energy (mechanical 
-|- gravitational) of a physical system is given by the ADM mass formula 

E ^ I ( - ,v^dl], (11) 



IQttG J V dx"- dxi 

where the integral is computed on a surface at spacelike infinity. We shall briefly review here the connection 
between this energy and the static gravitational potential. 

A generic static metric gp,i, can be written at spatial infinity in the form 

- - + ^ + '''' 

2M2G ^ M \ 

Performing the integral ( [Tl| ) with hij given by ( p^ ) one sees that M2 is the ADM energy ("total mass"). 
On the other hand. Mi is the mass observed by measuring the newtonian force at infinity. Substituting (^2|) - 
( p^ ) into Einstein's equations R^u = 0, it is easy to see that Mi — M2. This is a quite natural result ||l|. In 
other words it means that, according to special relativity conceptions, the source of the newtonian field is not 
only the mass, but also the energy density. For instance, in the gravitational collapse of a star a part of the 
gravitational energy is converted into kinetic energy and eventually this energy is employed to produce heavier 
elements from hydrogen or helium. If we disregard the radiation emitted into space, the newtonian field far 
away from the star remains unchanged during the whole process and the same holds for the ADM mass, which 
is a conserved quantity. 
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The gravitational potential energy can be found, by definition, assuming a static distribution of matter and 
computing the metric it generates at infinity. This has been done by Murchada and York 1^ for a spherical 
matter distribution of uniform density, using conformal transformations and a special formulation of the initial- 
value equations of General Relativity. For a sphere of (small) density p and unit radius they found the right 
newtonian gravitational binding energy, namely the ADM mass is given in this case by (reintroducing the radius 
R and the velocity of light c) 

Mtot = ^TiR'p - ^T^T^'p'Gi?^ + O(p'). (15) 

Remembering that A/tot also represents the effective source of the newtonian field, we see that the second 
term in (^|) gives rise to a deviation from the famous law which states the independence of the potential on 
the radius of the source. Nevertheless, this effect is usually unobservable, due to the very small factor c^^. 

It is also possible to find the following corrections to ([l5|), proportional to p'^, p*^, ... They denote the existence 
of general-relativistic corrections to the potential energy mim2G /r. For instance, the term proportional to 
would contribute to Mtot a term of the form 

AM oc ^p^G'^R' + Oip'^). (16) 

In the case of a source constituted by two pointlike bodies of masses mi and m2, kept at rest at a fixed 
distance r, the method of solution mentioned above is not applicable. From eq. ( |T5| ) we may infer that the 
ADM mass is given in this case by 

1 Gmim2 / 1 \ 

Mtot = mi + m2 tt h o ^ . (17) 

r \c 



There are, however, no relativistic corrections to the two-body potential coming from (^^, because the corre- 
sponding potential does not admit a continuum limit. For instance, if we try to integrate a potential of the form 
G^mJ^^mj^^/?'^ to obtain the term proportional to (? ^ we find that the binding energy of the sphere depends 
on the way it has actually been put together. This means that the correction proportional to comes from a 
three particles potential, and so on for p^ etc. 



3 Quantum formula for the potential energy on a fiat 
background. 

The same result we found in the preceding Section using the classical equations of motion can be obtained 
in a completely different way. It is known that the ground state energy of a system described by an action 
S'o['/>] — J d^x L{(p{x)) in the presence of external sources J{x) can in euclidean quantum field theory be expressed 
as 

h, fd[(j)]exp{-n-^[fd'^xL{(l){x))+fd'^X(j){x)J{x)]] 
E = lim log z 7 — 

T^oo T ^ ^ d[(t)\ei^p{^n-^ ^ d^xL{(j){x))} 

(18) 



= ^lim — ^ log ^exp ^ j d"^ x 4>{x) J {x)^ 

where, outside the interval {—\T, hT)^ the source has been switched off. 
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This equation has been proved exactly in perturbation theory Q for the case of a hnear focal coupling 
between the field and the external source, but it can be generalized assuming that in any case the vacuum-to- 
vacuum transition amplitude is given by 

^ / d[4^] exp { ~h-\So [4>] + Winter. J])} 

/d[(/.]exp{-n"i5oM} 

= (exp{n-l5i„ter.[</', J]}). (19) 

In fact, inserting a complete set of energy eigenstates we can write 

(0+|0-),7 = (0|e-^^/''|0} =^(0|e-^^/'"»(n|0) =^|(0|n}|2e-^"^/^ (20) 

n n 

The smallest energy eigenvalue En corresponds to the ground state, and in the limit T — > oo it dominates the 
sum. So taking the logarithm and multiplying by (—h/T) we obtain that energy. This is a well-known technique 
in QCD (see for instance 

In the case of a weak gravitational field quantized on a flat background, we may consider the source consti- 
tuted by two masses mi, TO2, placed at rest near the origin at a distance L each from the other (see eq. (^). 
The action of this system is 



S = - 

16 



-1712 



T 

^ J d'^Xy/^Rix) - mi y dh ^gf,^[x{h)]xt'{h)x''{h) 



dt2 ^gMt2)]r{t2)r{t2). (21) 



where the trajectories x^{ti) and 2/^(^2) of the particles with respect to the background are simply given by 

x^{ti) = (ti, 0, 0) ; y^(t2) - (t2, ^, 0, 0) . (22) 



2' ' / " ' ^' V 2 



So we have 

Ev!-a ~ hni x 



J d[h] exp {-h ^ S'Einst. + mi /J^: dti + hQo[x{ti)] + m2 J^t dt2 ^1 + /ioo[y(<2)] | 
log ^ ^ ' ' 



/ d[h] exp{-?i ^S'Einst.} 



^li^^_^log/exp -a-i 



mi / dti + hoa[x{ti)] +1712 dt2 + hoo[y{t2)] 



(23) 



We call this expression "B-G" (background-geometry) formula It has the property that the time T is 
referred to the background geometry. It follows in particular that the average is sensitive also to fields h which 
are pure gauge modes and carry no curvature. For further details on this point compare next Section and the 
discussion in Section 5. 

By standard perturbation techniques it is straightforward to see that for weak fields and to lowest order in 
G eq. (p3|) reduces to 



Eb-g = ™i + ^2 + lim^ log < 1 + I ^'^^^ I ^ {hoQ[x{ti)]hoo[y{t2)]) 
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1 'm\ui2 f ^ f ^ 
~ mi +1712+ lim - — — - — / dti / dt2 K^^^Jti ~ T2, L, 0, 0) 

T^OO 1 4 IT IT 

= mi+m2 (24) 

In the last step we have used the classical equation (|^). The next term in the perturbative series would be the 
first quantum correction, proportional to fiG^ . 

Eq. (|2^), like the corresponding ones in QED or QCD, has the physically appealing feature of showing how 
the force between the sources ultimately arises from the exchange of massless gravitons. However, let us make 
a closer comparison with electrodynamics. In that case the analogue of the functional integral which appears 
in the logarithm of ( ^3|) has the form 

^expjc, J dtiAo[x{ti)] - g J dt2Ao[y{t2)]^'^ . (25) 

(The two charges have been chosen to be opposite: qi = g, q2 = —g-) Reversing the direction of integration in 
the second integral and closing the contour at infinity, one is able to show that the quantity ( p5| ) coincides with 
the Wilson loop of a single charge thus giving a gauge invariant expression for the potential energy. 

In gravity this is not possible: we may imagine that an expression like ( p5| ) could be obtained in the first- 
order formalism (with ^0 replaced by the tetrad ej]), but the masses necessarily have the same sign, so the loop 
cannot be closed. 



4 General case. 

We aim now at giving a quantum formula for the potential which does not rely on a fixed background. We 
assume that a functional integral for euclidean gravity exists, denoted by 

Z = j d[g] exp {-ri-i^Einst.Lg]} (26) 

and require that all the field configurations in this functional integral are asymptotically fiat. 

In practice it is known that the euclidean Einstein action is not bounded from below, so one has to be careful 
in treating it. The action needs to be regularized in some way; in the actual lattice simulations, this happens 
thank to the term and to the measure (see and references). 

Since we have no background, the quantities T and L which enter in the definition of the source must now 
refer to the dynamical metric 5, configuration by configuration, in the functional integral (p6|). 

A first, simple recipe for generalizing the B-G equation ( p3| ) along this lines is the following. Let us consider 
two massive bodies (with mi = 7712, for simplicity) which are kept at a fixed invariant distance L from each 
other by some device and "fall" through a given field configuration g. Suppose that the motion is started and 
stopped in the asymptotically flat region. If T is the total proper time measured along the trajectory of the 
first mass (with T ^ L), the proper time measured along the trajectory of the other mass will be, say, T + a[g]. 
Then the action of this system leads, after integration over all field configurations, to the energy 

Et-g - lim -^\og{exp{~h-^mi2T + a[g])}) 

= 2m + ^lim -^log(exp{-?i"^TOa[5]}) . (27) 

We call this expression "T-G" (total geometry) formula. Unlike the B-G formula (eq. ([23|)), it refers only to 
geometrical, invariant quantities. 
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If we suppose that the field configurations contributmg to the average are not too smgular, then a is small 
and we may expand the exponential, finding (note that (a[g]) obviously vanishes by symmetry) 

ET^G = 2m+\im^~ — {a'[g]) + ... (28) 

Thus the interaction energy is always negative in this case. Also, eq. ( p8|) gives an intuitive picture of the 
gravitational attraction: it arises because keeping two masses at some fixed distance increases their action, due 
to the vacuum fluctuations of the gravitational field. 

Now, we turn to a more careful description of the whole procedure above, and eliminate the (apparent! - 
thank to the average) asymmetry between the two bodies, assigning the total proper time T to the geodesic of 
their center of mass. 

Let us suppose that a field configuration is given. We consider a geodesic line of length T, which starts at 
an arbitrary point in the "past" asymptotically flat region with unit timelike velocity. 
To fix the ideas, this curve could be written in its first part as 

e{r)^ (^-\t + t, 0,0, 0<r<f, (29) 

where r is the proper time measured along the curve and we have chosen the remaining coordinates of the 
starting point to be equal to (0, 0, 0) (this is an irrelevant arbitrariness, since at the end we shall integrate over 
all the configurations of the field). As usual, T denotes a very long time interval. After a time ~ f the curve 
enters the region of spacetime where the gravitational field is non vanishing. It continues as a geodesic, which 
means that £,^{t) satisfies the equation 

KMrW{r)t{T)+e{T)=0, (30) 

where FJ^^, is the Christoffel symbol of the metric. The curve terminates at t = ^T, again in the flat region. 

Let us then take in the initial point ^'^(0) a unit vector <7^(0), orthogonal to ^'^(0) (for instance, in our 
example, q^{0) = (0, 1, 0, 0)), and define a vector q^{T) along the curve ^^{t) by parallel transport of q^{0). 
We remind that ^'^{t), being the tangent vector of a geodesic, is parallel transported along the geodesic itself, 
and that the parallel transport preserves the norms and the scalar products. Then the following relations hold 
along the curve 

e{r)C{T)g,Aar)] = -1; (31) 
q'{r)q''{r)9,Aar)] = 1; (32) 
eW{T)g,Aar)] = 0. (33) 

Next we consider two masses mi, m2, and a length L which we may regard as infinitesimal, compared to 
the scale T. We assume that the two masses follow the trajectories x^{t) and y^ir), respectively, given by 

xf^^r) = enr)-Liq^(r); (34) 
y^ir) = e'{r) + L2qnr), (35) 

where Li and L2 are two positive lengths such that 

Li + L2 — L and — wiii + 7712^2 — 0. (36) 
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The physical meaning of the preceding geometrical construction is apparent: it represents an observer which 
falls freely in the center of mass of the system composed by mi and m2, while holding the two masses at rest at 
a distance L each from the other. This is the generalization of the source introduced in eq. (|2^) that is naturally 
dictated by the equivalence principle. 

We notice that if the two masses were allowed to fall freely in the field, they would not keep at a constant 
distance from each other. In fact, as it is well known from the so-called geodesic deviation equation, the distance 
between two neighboring geodesies varies according to the sign of the curvature in the region they are traversing. 

We can reparameterize the two curves x'^(r) and y''(T) introducing their proper times ri and T2, respectively. 
The ratio between the proper time ti and the proper time r is given by the equation 

dn ^ -gA<Ti)\i^{'ri)i''{Ti)dT. (37) 

An analogous relation holds for T2. We agree to adjust the fmiction ti(t) in such a way that ti(0) = 0. 
Then we shall denote Ti(-iT) = -\T[, ti{\T) = \T\' and Ti = \{Tl + For flat geometries we have 

T[=T'^ = Ti=T. Analogous relations hold for T2. 
The energy can then be written as 

Et^g = ^lini^ log (exp {~h~\miTi + maTa]}) . (38) 



In order to make contact with eq. (|2^), suppose now to take mi = m2 and to "attach" smoothly the ends 
of the curves x{t) and ?/(t) to the geodesic ^(r). The corresponding change of length will have no effect on eq. 
(^), because it is in no way a function of T . Since a geodesic has minimal length with respect to neighboring 
curves, we have, for the total times Ti and T2 as functions of L 

^ = l + la[g]L^ + \b[g]L^ + ... (39) 

Y = l + \a[9]L^-\b[g]L^ + ... (40) 

where a[g\ > 0, while b has no definite sign. From this we can find the difference between Ti and T2, which is 
the analogue of the quantity a[g] of eq. (p7|): 

g = l + + that is ^^b[g]L\ (41) 



5 Discussion. 

The B-G formula ( p3| ) and the more general T-G formula (|^ (which is the one actually implemented in the 
lattice simulations) must be consistent in the limit of small G. In this Section we are going to clarify this point 
and work out some consequences. 

We first recall that expanding the euclidean functional integral ( p6| ) from fiat space in perturbation theory, 
one finds that to order HG all field configurations are pure gauge modes carrying no curvature. In fact, the 
Wilson loops of the connection vanish to this order, and this cannot be an average effect, as they have the form 
W — {9i + 62}, where 61 and 62 are the two angles which describe the 0(4) rotation of vectors by parallel 
transport around the loops (see [Tit). 

""^This also suggests that gravitons are, physically, quite fictitiuos objects. 
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As we pointed out in Section 3, the B-G formula is sensitive to such pure gauge modes, and gives in fact to 
perturbative order hG the right result for the potential. On the other hand, the T-G formula would clearly give 
no potential energy in flat space: if there is no curvature, the difference between the lengths of the trajectories 
of the two masses vanishes. So the two equations "work" in a very different way. 

In lattice theory, a version of the T-G formula is employed. In their non-perturbative simulations, Hamber 
and Williams Q trace two parallel lines with reference to a (dynamical) Regge lattice and then compare the 
lengths li of the two lines. As we mentioned in the Introduction, the first results are positive. 

In general, the lattice formulation seems to be much more appropriate for handling quantum gravity than the 
perturbative expansion. It also has the property of being automatically coordinate invariant, and of generating 
dynamically the flat space limit. The perturbation theory on flat background makes sense instead as an effective 
theory for small G, where G is the effective value of the gravitational constant as measured from the macroscopic 
Newton force. 

In the phase space of lattice theory, specified by the parameters a (the coupling to the i?^-term) and GBare, 
there is a continuous line which separates a physical "smooth" phase and a collapsed "rough" phase |^ . The true 
continuum theory is the limit of lattice theory on a point of this line, chosen in such a way that the effective 
long-distance Newton constant, computed through the lattice version of T-G formula, equals the measured 
value. 

Then we may argue that a "correspondence principle" between the B-G and T-G formulas requires that 
Gsarc is not as small as the effective value. Namely, for the region of very small Gearc first order perturba- 
tive computations are reliable, so we know there is no curvature in that regime, and consequently the lattice 
evaluation of T-G formula will give G — exactly (or better, apart from terms of order Gq^j.^). On the con- 
trary, if GBarc is not so small, its effective value can emerge from lattice theory with no contradiction with 
the background computations. Whether this is really what happens, it should become clear as the simulations 
proceed. 

It is a pleasure for me to thank D. Maison for the warm hospitality in Munich, and the A. Von Humboldt 
Foundation for financial support. Also part of this work was made at M.I.T., Center for Theoretical Physics, 
and the help of R. Jackiw is gratefully acknowledged. 
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